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ABSTRACT 

We present a new approach for calculating tunnelling amplitudes from a 
non-localised initial state. Generalising the matching conditions and equations 
of motion to allow for ‘complex’ momentum permits a description of tunnelling 
in the presence of so-called classical motion. We comment on possible applica- 
tions of the method. 
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That quantum mechanical probability density can permeate into regions in which a 
classical particle is forbidden was (and perhaps still is) one of the more intriguing features of 
the Schrodinger equation. Rapidly however, it was turned into one of the major successes of 
the early wave mechanics in providing a successful theoretical description of alpha particle 
emission in radioactive decay. Since then, tunnelling has found many applications, yet 
our understanding of many aspects of the phenomenon remains limited. Motivated by 
difficulties 1 in calculating tunnelling amplitudes in recent two-field inflationary models 2- " 4 , 
we re-examined the quasi-classical approximation to the Schrodinger equation to see if 
one could incorporate more complicated evolution of the wave function without losing the 
simplicity of that approach. We found that the problem reduces to two key issues: that 
of incorporating complex momentum into a quasi-classical description, and the related 
problem of matching conditions between different asymptotic regimes. Here we will show 
how to resolve these questions, illustrating some results from a test potential. A more 
detailed exposition may be found in a companion paper 5 . 

First let us summarise the current status for calculating tunnelling amplitudes. Most 
tunnelling calculations proceed via the quasi-classical approximation 6 , which extracts the 
leading order “classical” behaviour of the system. If we assume the wave function takes 
the form 

*(x,<) = e i<r ^/ h e~ iEt / h , (1) 

then the Schrodinger equation gives 


(Vtr) 2 - — VV = E - tf(x) 
2m' 2m 


( 2 ) 


as the equation of motion for a . The quasi-classical approximation drops the 0{h) term 
in (2), which clearly requires that |V<r| 2 » fi|V 2 <r|. 

In tunnelling, we need to solve the Schrodinger equation in regions where U > E, 
therefore we set a = ig and obtain 


(3a) 


= 4 > 


(Vs) 2 = 2 m(U - E) 

g= f ' y/lMV - E)dl, 


(3b) 



where 1 is the path which minimises the integral and X/ is the point of emergence of that 
path into the classical regime. This geometric solution to the problem in more than one 
dimension was first obtained by Banks, Bender and Wu 7 for the case of tunnelling from 
a localised state, they called 1 the escape path . The calculation of 1 and g in practise 
is facilitated by the Euclidean time prescription whereby we associate Vg with p e , a 
Euclidean momentum, and (3a) becomes a H ami ltonian problem of particle motion in the 
(inverted) potential -17. It is then a classical dynamics problem to calculate the trajectory 
x(r), which interpolates between the initial position and the position of emergence from 
the barrier. 

This approach is powerful, however it does have one crucial restriction, namely, it 
requires that space be divided into regions in which V<r is entirely real or entirely imaginary 
so that we can identify p = Vcr with some x(“time w ). Clearly we then need x to vanish on 
the boundaries between these regions. In the context of the escape path, we need x 2 — ♦ 0 
at each end of the path. In general it requires two initial conditions to fix the start of 
the escape path, which uses up the requisite number of boundary conditions, leaving no 
additional freedom for the end point, which in general will not satisfy x 2 = 0. Another 
related problem is that we have no allowance for the transport of real momentum under a 
barrier, which we know definitely does occur in the case of a continuous symmetry. There 
is also the problem of picking an initial position from which to integrate, since we can only 
localise a paxticle at the expense of information about its momentum. We clearly need a 
more general picture of tunnelling. 

In order to solve the Schrodinger equation when the momentum is complex, we rewrite 
(2) in terms of the real variables / and g } where a = / + ig: 

(V/) 2 - (■ V<?) 2 + hV 2 g = 2 m(E - U) (4a) 

2 V/.Vtj - hV 2 f = 0. (46) 

In the quasi-classical approximation we neglect the 0(h) terms in these equations, this will 
be valid unless 

(V/) 2 ~ hV 2 f 

( 5 ) 

(V< 7 ) 2 ~ hV 2 g 
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i.e. unless J V/[, |V^| is small or |V 2 /|, |V 2 ^| large, this latter situation occurs at caustics 
of the motion. The simplest way to visualise this is to consider a plane wave scattering 
off some barrier. The integral curves of V/ trace out the path that a family of classical 
particles would follow in that potential. Clearly, neighbouring, initially parallel, trajecto- 
ries will at some stage cross, and at such a place, referred to as a caustic, V 2 / becomes 
unbounded. 

Having identified the boundary between different asymptotic regions, we now focus on 
the turning surface to obtain the matching conditions. We first consider a local coordinate 
system along the caustic {x,i/*} where x is orthogonal distance from the caustic, and the 
y % are cartesian coordinates in the surface of the caustic. Provided the caustic is not too 
strongly curved, the potential is roughly a function of x, and we may deduce that the true 
(continuous) solution of the Schrodinger equation is 

( 6 ) 

where Ai(x) is the Airy function. That the curvature of the caustic not be too large trans- 
lates to a bound on V\U: noting that (6) implies that the quasi-classicai approximation 

breaks down when |x| < h 2 t*(2mU, z J"” 1 / 3 , we obtain a limit of V\U <C . Equa- 

tion (6) gives us the inital conditions for / and g for integrating out beyond the caustic. 
Clearly the parallel momenta, p. are conserved across the boundary, i.e. V t /, con- 
served. The orthogonal components are determined on either side of the boundary by using 
equation (4a) and demanding that real (V/) and imaginary (Vj) momenta be orthogonal 
(equation(4b)). This fixes the initial conditions for the momentum on the other side of the 
boundary. We now need to propagate / and g out from the boundary, i.e. solve (4a, b). 

In order to facilitate the solution of (4a) we assume E If, which implies that V <7 is 
dominant under the barrier, however, we do not wish to neglect V/ as an O(fe) correction, 
otherwise we would use existing techniques. Instead we want to consider a situation where 
p <1, and iteratively solve 

(V 5 ) 2 = 2 mU - (2 mE - (V/) 2 ) 


V>(x,y‘) - e* p ' yi f h Ai h^ 2/ \2mU yz )*/ 3 


x + 


E - 17(0) - 

u7~ 


2m 


(7a) 


bearing in mind that (V/) 2 , E are of the same order, and small compared with (Vp) 2 , U. 

The first step is to find the leading behaviour, that is, to solve (Vy) 2 = 2mI7. In 
order to do this, we use existing techniques. We regard the solution under the barrier as 
being composed of a family of escape paths {l (a) } originating from the turning surface at 
(y(<»)’®)* propagating under the barrier until a caustic is reached. Along each escape 
path we have a solution for g determined by (3b) with E = 0. This will give us a leading 
order solution to the problem. Now we use the information from (7b), that / is constant 
along integral curves of g , to set / = /(<P( Q )( X )) where ¥>( a )( x ) = !/(„) 8X6 the equations 
for the family of escape paths l( a ). Finally, we input this solution for / back into (7a) to 
obtain the correct form of g to order EjU. Once we have the solution under the barrier 
we may match across into the far ‘classical’ regime and follow a similar procedure to the 
above to complete the solution of the Schrodinger equation. 

To recap: this method solves the Schrodinger equation in the stationary quasi-classical 
approximation for tunnelling in potentials with non-localised minima. The restrictions on 
the type of potentials it can deal with are that E U under the barrier and that the 
transverse derivatives of U at the matching boundary be bounded by h ^ . This 

method is therefore ideal for problems involving scattering off a ‘wall’. 

As a simple illustration of the method, consider a plane wave, e , * c - x /^ ) scattering off 
a wedge potential in two dimensions, C7(x) = V0(i)0(a — x — ey), see fig. 1. Here, our 
initial conditions at x = 0 are g = 0, / = 5>* = \/2mF and /,* — 0- Integrating 

out is straightforward since the potential is constant and gives g = y2mV x. Imposing 
/ =const. along integral curves of <7, which in this case are y = y., simply gives / = k } y, 
then inputting this back into (7a) gives g = y/2mV — k*x = kx. In this example our 
iterative procedure under the barrier is now complete. At the far edge of the barrier the 
appropriate boundary conditions are now non-standard, the momentum parallel to the far 
edge being complex: p |( = (p, — te«)T, where T is the tangent vector to the far edge 
of the barrier. Orthogonality of the real and imaginary momenta and (7a) requires that 
px = (p t +ieK^)N+0(e 2 ). Setting 


( = x - a + ey , r 1 = y - c(x - a) 


( 8 ) 
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(see fig. 2) we may integrate out from the boundary to obtain 


. p, , m 

J.„ =Ka- cktj H 

Pi 

Note some interesting features about this solution. The incoming and outgoing real mo- 
menta axe not parallel, although they subtend the same angle with each side of the barrier. 
Note too that g is spatially dependent, decreasing in regions where the barrier is thinner. 
However, the really curious features of scattering off this potential axe most clearly ex- 
hibited if we consider a Gaussian wave packet scattering off the barrier. For the sake of 
definiteness we take the momentum to be peaked at (Pi,p t ) with spread A ft. 2 , hitting the 
origin at t = 0. Such an incident wave packet might reasonably be supposed to represent 
a particle scattering off the barrier. For the trajectory of the outgoing peak we obtain 


p,(l — 2aAh//c) 2A h 


+ 1 + S(1 ~2a\h/ K y 


with a transmission amplitude of e~** out /*. In other words, a particle ‘incident’ on the 
barrier atf = x = y = 0 emerges at = - ^(1 + (1 — 2a\h/ k)*). One can also calculate 

the time of emergence, t c = (1 — 2aAft//c) 3 . The particle emerges considerably 

‘downstream’ from where one might expect it, neither at x = (a, 0), nor at x = (a, ea), but 
at x = (a—erjc y rf e ). Furthermore, the time at which it emerges is potentially large compared 
with h and even, if p 7 < 0, potentially negative! In fact, both these peculiarities arise 
from a simple physical reason, the fact that the most energetically favourable time for the 
Gaussian to tunnel is not necessarily when the peak hits the barrier. Tunnelling amplitudes 
depend exponentially on the size of the barrier, therefore it is more favourable to tunnel 
where the barrier is thinner, hence the dependence of the transition amplitude on x out . On 
the other hand, the probability density along z = 0 is damped by an exponential factor 
depending upon how far away that point is from the peak of the Gaussian. Cleaxly there 
will be a pay off between these two factors which may mean that it is more energetically 
favourable for the fringe of the Gaussian to tunnel, rather than its peak. 
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Apart from giving an initial step in the problem of calculating two-field tunnelling 
amplitudes, our method (and illustrative example) may also have applications in other 
areas of tunnelling theory. For instance, one of the outstanding problems in tunnelling 
is defining a unanimously acceptable tunnelling time. From the uncertainty principle, we 
know such a time, tx , should be of order h , but when using stationary states in a time 
independent potential, finding tx is not so straightforward, indeed quite a controversy 
exists®. Our example suggested that one of these definitions, the extrapolated phase time 
is not a good definition of tunnelling time as it stands. Unfortunately, we can shed no light 
on what is! 

While we have only applied this method to a s ma l l range of examples 5 , the results 
we have so far obtained are encouraging: the modifications seem to give quite different 
qualitative pictures for scattering. The esssential step forward from the Banks, Bender 
and Wu approach was to identify a general boundary between different WKB regimes. It 
can only be hoped that a similar generalisation of Coleman’s field theoretic procedure 1 is 
possible. 
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